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1 INTRODUCTION 



ABSTRACT 

Numerical simulations of the magnetorotational instability (MRI) with zero initial net flux in 
a non-stratified isothermal cubic domain are used to demonstrate the importance of magnetic 
boundary conditions. In fully periodic systems the level of turbulence generated by the MRI 
strongly decreases as the magnetic Prandtl number (Pm), which is the ratio of kinematic 
viscosity and magnetic diffusion, is decreased. No MRI or dynamo action below Pm = 1 is 
found, agreeing with earlier investigations. Using vertical field conditions, which allow the 
generation of a net toroidal flux and magnetic helicity fluxes out of the system, the MRI is 
found to be excited in the range 0.1 < Pm < 10, and that the saturation level is independent 
of Pm. In the vertical field runs strong mean-field dynamo develops and helps to sustain the 
MRI. 
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The realization of t he astrophysical signific ance of the magnetoro- 
tational instability ("Balbus & Hawlev 1991), first discovered in the 
context of Couette flow (Velikhov 1959; Chandrasekhar 
[ seemed to resolve the long-standing problem of the mechanism 
• driving turbulence in accretion disks. Early numerical simulations 
produced sustained turbulence, large-scal e magnetic fields and out- 
, ward angular momentum transport (e.g. Brandenb urg et ^ 1 19951 : 
■ lHawlevetal1 IT995). These results also showed that a significant 
qualitative difference exists between models where an imposed uni- 
form magnetic field is present as opposed to the situations where 
such field is absent: the saturation level of turbulence and angular 
momentum transport are substantial ly higher when a non-zero ver- 
tical net flux is present (e.g. Brandenburg et alll 19951 : IStone et aEI 
1 19961) . Also the presenc e of an imposed n et toroidal field seemed 
to enhance the transport JStoneet all 19961) . 

In the meantime, a lot of numerical work has been done 
with zero net flux setups that omit stratification and adopt fully 
periodic or perfectly conducting boundaries in order to study 
the saturatio n behaviour of the MRI in the simplest possible 
setting (e.g. Fromang & Papaloizou '2007|: iFromang et al.l 1 20071 : 
[Liljestrom et al. 2009; Korpi et al. 2010). Due to the boundary con- 
ditions, the initial net flux in conserved and no magnetic helicity 
fluxes out of the system are allowed. The results of these investi- 
gations have shown that as the numerical resolution of the simula- 
tions increases, or equivalently as the explicit diffusion decreases, 
the level of turbulence and angular momentum transport transport 
decrease, constituting a convergence problem for zero net flux MRI 
jpromang et al.ll2007l) . Runs with explicit diffusion show that sus- 
taining turbulence becomes increasingly difficult as the magnetic 



Prandtl number, Pm = v/ri, wh ere v is the viscosit y and ri the 
magnetic diffusivity, is decreased (Fromang et al. 2007). Currently 
the convergence problem is without a definite solution. It has been 
suggested that this issue could be related to the Pm-dep endence 
of the fluctuation dynamo (e.g. ISchekochihin et al]|2007l) . It has 
even been argued that the MRI in periodic zero net flux systems 
would vanish in the limit of large Reynolds numbers and that a 
large-sca le dynamo wou ld be needed to sustain the MRI and tur- 
bulence dVishniadlioogl) . Notably, large-scale dynamos have no 
problems operating at low magnetic Prandtl numbers as long as 
the relevant Reynol ds and dynamo numbers exceed critical values 
( lBrandenbur3l2009l) . 

From the point of view of mean-field dynamo theory 
terandenburg & Subramanianll2005l) , systems with fully periodic 
or perfectly conducting boundaries are rather special. In such 
closed setups magnetic helicity, defined as a volume integral of 
A ■ B, where A is the vector potential and i? = V x A is the 
magnetic field, is a conserved quantity in ideal MHD. In the pres- 
ence of magnetic diffusion, magnetic helicity can change only on a 
timescale based on microscopic diffusivity, which is usually a very 
long in any astrophysical setting. Such a behaviour, which has been 
captured in numerical simulations (^.randenburg 2001), is well de- 
scribed by simple mean-field models taking into account ma; 
netic helicity conservation (e.g. iBlackman & Brandenburg! 1200: 
This would mean that generating appreciable large-scale mag- 
netic fields, which are possibly vital for sustaining the MRI, 
can take a very long time. Furthermore, the saturation value 
of the mean magnetic field decreases inversely proportional to 
the magnetic Reyn olds number (e.g. ICattaneo & HughesI 1 19961 ; 
iBrandenburd I2001I) . In dynamo theory this detrimental effect to 
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the large-scale dynamo is kn own as the catastrophic quenching 
jVainshtein & Cattaneoll 19921) . 

The situation, however, changes dramatically if magnetic 
helicity flux ou t of t he system is allowed. In particular, the 
IVishniac & Choi bOOlh flux, which requires large-scale velocity 
shear to be present and flows along the isocontours of shear, is a 
potential mechanism that can drive a magnetic helicity flux out of 
the system and alleviate catastrophic quenching. Indirect evidence 
for its imp ortance exists from convec tion simulations in a shearing 
box setup jKapvla et alj|2008l , l2010bi) . where dynamo excitation is 
easier in systems with boundaries that allow a net magnetic helicity 
flux. However, these results can be explained by a somewhat higher 
critica l dynamo number in the perfect conductor case ( Kapvla et al. 
l20I0bl) . which is a purely kinematic effect. More dramatic differ- 
ences between different boundary conditions are seen in the non- 
linear saturation regime, with strong que nching of large-scale mag- 
netic fields in the perfect conductor case ( Kapv la et al .l2010 b!). The 
reason for this behav iour is not yet clear , espec ially in light of re- 
cent results of IHubbard & Brandenburd who failed to find 
evidence of the Vishniac-Cho flux in a numerical setup similar to 
ours. 

In the present paper we demonstrate that the boundary condi- 
tions play a crucial role for the excitation of the MRI and the associ- 
ated large-scale dynamo. Following previous work that has shown 
that open boundary conditi ons allow more efficient dynamo action 
jKiiDvla e t al. 2008, 20I0bl) . we model a system that is isothermal, 
non-stratified, and the magnetic field has a zero net flux initially. 
We then apply vertical field boundary conditions which allow a 
magnetic helicity flux through the vertical boundaries by letting 
the magnetic field cross them. We show that if the MRI is excited, 
a large-scale dynamo is also excited and that the saturation level 
of the turbulence, large-scale magnetic field, and angular momen- 
tum transport are essentially independent of Pm. This is contrasted 
by periodic simulations where we find a strong Pm-dependence in 
accordance with earlier studies. Our results also suggest that for a 
given Pm the results (level of turbulence and angular momentun 
tran sport) are indep endent of the magnetic Reynolds number (see 
also lFroman 3l2010h . 

The remainder of the paper is organised as follows: in Sect.|2] 
we describe our model, and in Sect.[3]and|4l we present our results 
and conclusions. 



2 THE MODEL 

In an effort to keep the system as simple as possible, we assume 
that the fluid is non-stratified and isothermal. The diffusion pro- 
cesses are modeled with explicit Laplacian diffusi on operators with 
consta nt c oefficients. A simila r model was used bv lLiliestrom et all 
(12009') and'Ko rpietal] ( l2010l) . although in these models higher or- 
der hyperdiffusive operators were used instead of the Laplacian 
ones. The computational domain is a cube with volume H'^ — 
(27r)^. We solve the usual set of hydromagnetic equations in this 
geometry 

^ = -SAyX-{VU)^A~riiioJ, (1) 
V\np 



^ - -V-LT, 
Vt 

VU 2 



(2) 



where V/Vt — d/dt + {U + u''^'^) ■ V is the advective time 
derivative, A is the magnetic vector potential, i? = V x A is the 
magnetic field, and J = fi^^'V x B is the current density, fio is 
the vacuum permeability, rj and v are the magnetic diffusivity and 
kinematic viscosity, respectively, p is the density, U is the velocity, 
and ft — Qo{0, 0, 1) is the rotation vector. The large-scale shear is 

given by U"'"' — (0, Sx, 0), with q = —S/Qo = 1.5, correspond- 
ing to Keplerian rotation, in all runs. We use isothermal equation 
of state p = c^p, characterised by a constant speed of sound, Cs. 
In the present models we choose the sound speed so that the Mach 
number remains of the order of 0. 1 or smaller in order to minimize 
the effects of compressibility. The rate of strain tensor S is given 
by 



(4) 



where the commas denote spatial derivatives. The initial magnetic 
field can be written in terms of the vector potential as 



A — Ao cos(fcA3;) cos{kAz)ey, 



(5) 



where the amplitude of the resulting magnetic field that contains x 
and ^-components is given by Bq = A;a^o. We use kA/ki = 1, 
Qo = I • lO'^Csfci, and ylo = ^ ■ ^ fj-oPoCsk^^ in all models. 

The values of /ca, and Ao are selected so that both the 
wavenumber with the largest growth rate, fcmax ~ Q,o/ua = 2, 
where ua ~ BoZ-^MoPo is the Alfven velocity, and the largest 
unstable wavenumber, fccrit ~ ^/^qkmax ~ 3.5, are well resolved 
by the grid. The other condition for the onset of MRI, namely /? > 
1, where /3 — 2pqp/Bq is the ratio of thermal to magnetic pressure, 
is also satisfied as /3 = 1800 for the maximum values of the initial 
magnetic field. 

We use the PENCIL COD^ which is a high-order exphcit fi- 
nite difference method for solving the equations of compressible 
magnetohydrodynamics. Resolutions of up to 512^ are used, see 
Figure[T]for a snapshot of a high resolution run. 

2.1 Boundary conditions 

In all models the y-direction is periodic and shearing- 
periodic boundary condit ions are used for the s-direction 
dWisdom & Tremainelll988l) . On the ^-boundaries we use two sets 
of conditions. Firstly, we apply periodic boundaries (denoted as 
PER). 

Secondly, we apply a vertical field (VF) condition for the mag- 
netic field, which is fulfilled when 



Bx = By 



B. 



0, 



(6) 



at the z-boundaries. In this case we use impenetrable, stress-free 
conditions for the velocity according to 



(7) 



The novel property of the VF conditions is that they allow a net 
toroidal flux to develop and allow magnetic helicity fluxes out of 
the domain. 



2.2 Units, nondimensional quantities, and parameters 

Dimensionless quantities are obtained by setting 

fci = Cs = po = Mo = 1 , (8) 



+ -{J X B + -V ■2upS), (3) 
P 



^ http ://pencil-code . googlecode .com 
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Table 1. Summary of the rans. The Mach number (Ma) is given by equation )14t . -Brms = Brms/^?cq, and Bi = \J 3^ /Bcq, where Bcq is defined via 
equation )13t . Kxy = Rxy / {^oH)'^ and M^y = {potJ-o)~^ Mxy / {^oH)'^ , where R^y and Mxy are computed from equations {TTJ and U8t . respectively. 
Finally, ags is given by equation M(>\ . 



Run 


grid 


Cm 


Rm 


Pm 


Ma 






R 


R [10 — 31 


M flO— 31 

^^^xyY^^ J 




BC 


AO 


128^ 


5 ■ 103 




5 




- 


- 


- 


- 


- 


- 


PER 


Al 


128^ 


10-* 


208 


5 


0.021 


2.09 


0.09 


0.51 


0.315 


-2.162 


2.477 ±0.270 


PER 


A2 


128^ 


1.5 ■ 10"* 


326 


5 


0.022 


2.04 


0.08 


0.54 


0.378 


-2.337 


2.715 ±0.208 


PER 


A3 


256^ 


3 ■ lO'' 


706 


5 


0.024 


1.92 


0.07 


0.35 


0.389 


-2.564 


2.953 ±0.338 


PER 


A4 


256^ 


3 ■ 10"' 


377 


2 


0.013 


1.78 


0.04 


0.31 


0.102 


-0.626 


0.728 ± 0.212 


PER 


A5 


256^ 


6 ■ 10" 


625 


2 


0.010 


1.83 


0.04 


0.33 


0.079 


-0.441 


0.520 ±0.074 


PER 


A6 


256^ 


3 ■ 10'' 


211 


1 


0.007 


1.28 


0.02 


0.34 


0.011 


-0.075 


0.086 ±0.022 


PER 


A7 


256^ 


6 ■ 10'' 


348 


1 


0.006 


1.57 


0.02 


0.31 


0.015 


-0.088 


0.103 ±0.014 


PER 


BO 


128^ 


1.5 ■ 10" 




20 
















VF 


Bl 


128^ 


1.5 ■ 10" 


557 


10 


0.037 


2.76 


0.12 


2.30 


0.866 


-4.726 


5.592 ±0.325 


VF 


B2 


128^ 


1.5 ■ 10* 


530 


5 


0.035 


2.06 


0.12 


1.18 


0.899 


-4.802 


5.702 ±0.299 


VF 


B3 


128=* 


1.5 ■ 10* 


632 


2 


0.042 


2.33 


0.12 


1.91 


1.140 


-4.577 


5.717 ±0.071 


VF 


B4 


128^ 


6.0 ■ 103 


307 


1 


0.051 


1.95 


0.13 


1.54 


1.519 


-5.463 


6.982 ± 0.909 


VF 


B5 


128^ 


1.5 ■ 10* 


637 


1 


0.042 


2.24 


0.12 


1.82 


1.164 


-4.422 


5.586 ±0.526 


VF 


B6 


2563 


3.0 • 10* 


1242 


1 


0.041 


1.77 


0.11 


0.97 


1.018 


-5.094 


6.111 ±0.560 


VF 


B7 


2563 


1.5 ■ 10* 


687 


0.5 


0.046 


1.69 


0.12 


1.04 


1.154 


-4.988 


6.142 ±0.636 


VF 


B8 


512^ 


1.5 ■ 10* 


719 


0.2 


0.048 


1.55 


0.11 


0.87 


1.111 


-5.076 


6.187 ± 1.068 


VF 


B9 


5123 


1.5 ■ 10* 


897 


0.1 


0.060 


1.78 


0.12 


1.39 


1.680 


-6.148 


7.828 ± 1.335 


VF 



where po is the mean density. The units of length, time, velocity, 
density, and magnetic field are then 



[x] = fcr' , M = (csfci)-^ , [u] = 



(9) 



po , [n\ = y /ioPoCs • 

The simulations are controlled by the following dimensionless pa- 
rameters: the magnetic diffusion in comparison to viscosity is mea- 
sured by the magnetic Prandtl number 

Pm=-. (10) 

The effects of viscosity and magnetic diffusion are quantified re- 
spectively by the parameters 

Cm - — 

rik\ ' Pm uk-^ 

We also define the fluid and magnetic Reynolds numbers 



(11) 



Re = 



i/ki ' 



Rm = 



rjki 



= PmRe, 



(12) 



where Mrma is the root-mean-square (rms) value of the velocity, 
better decribing the nonlinear outcome of the simulations. Further- 
more, we often measure the magnetic field in terms of the equipar- 
tition field which is defined via 



-Boq = \/ /io(pU?ms), 



(13) 



where the brackets denote volume averaging. A convenient mea- 
sure of the turbulent velocity is the Mach number 

^rms 



We define the mean quantites as horizontal averages 



■-Lx/2 



Ly/2 



Fi{x,y,z,t)dxdy. 



(14) 



(15) 



-Ly/2 
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Figure 1. Velocity component Ux from the periphery of the domain for 
Run B9 with Pm = 0.1, Cm = 1.5 • 10*, and Re 9 • 10^. See also 
http : / / www .helsinki . fi/ '-^kapyla/movies . html for anima- 
tions. 



Often an additional time average over the statically saturated state 
is also taken. The size of error bars is estimated by dividing the 
time series into three equally long parts. The largest deviation of 
the average for each of the three parts from that over the full time 
series is taken to represent the error. 
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Figure 2. Mach number defined via equation U4t for Runs A3-A7. The 
thick solid line shows the Mach number for Run B9 with Pm = 0.1 and 
VF boundaries. 
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Figure 3. Mach number (upper panel) and magnetic energy (lower panel) as 
functions of magnetic Prandtl number for periodic (triangles) and vertical 
field (diamonds) boundary conditions. The magnetic field is normalised by 
the rms value of the initial field. 



3 RESULTS 

We perform two sets of simulations listed in Table[T]where we use 
either periodic (Set A) or vertical field (Set B) boundary conditions. 
In Set A, Runs A0-A3 were started with the initial conditions de- 
scribed in Section [2] whereas Runs A4-A6 were continued from a 
snapshot of Run A3 in the saturated state, see Fig. |2l Run A7 was 



continued from a snapshot of Run A6 with a two times lower dif- 
fusivities at roughly 150Torb where Torb = 27r/f2o is the orbital 
period. The minimum duration of the runs in Set A is lOOTorb- 
Runs in Set B were all started from scratch and typically ran a sig- 
nificantly shorter time than those in Set A, e.g. ~ SOTorb in the 
low-Pm cases (see Fig. |2j, because final saturation occurs much 
faster. 



3.1 Saturation level of the MRI 

3.1.1 Periodic case 

Earlier studies have shown that exciting the MRI in a periodic 
zero net flux system becomes increasi ngly harder as the magnetic 
Prandtl number is decreased jFromang et al. 2007.) . Furthermore, 
the saturation level of turbulence has been reported to decrease as a 
function of Pm. This has been conjectured to be associated with the 
difficulti es of exciting a small-sca le or fluctuation dynamo at low 
Pm (e.g. lSchekochihin et al.l2007l) . It is, however, unclear how the 
saturation level of the small-scale dynamo is affected by this. It is 
conceivable that at magnetic Reynolds numbers close to marginal it 
takes a long time to reach saturation and that the current simulations 
have not been run long enough. On the other hand, if catastrophic 
quenching is t o blame, the mean magnetic field should decrease 
as Rin"^ (e.g. lBrandenburg & Subramanianll2005i . and r eferences 
therein ). A further possibility is the scenario suggested bv lVishnia3 
( l2009h : in the absence of an outer scale for the magnetic field, the 
microscopic diffusivities determine the minimum lenght scale of 
MRI, which leads to turbulence intensity decreasing proportional 
toRm-2/^ 

We study this issue by performing runs keeping Pm fixed and 
increasing the Reynolds numbers. We find that the saturation level 
of turbulence, measured by the Mach number and root mean square 
value of magnetic field, are unaffected when Cm is increased by a 
factor of three for the case Pm = 5 (Runs A1-A3) and by a fac- 
tor of two for the cases Pm — 2 (Runs A4-A5) and Pm — 1 
(Runs A6-A7), see Table \T\ and Fig. |3] Furthermore, the Mach 
number and rms magnetic field, normalised with the rms value 
of the initial field, increase roughly linearly with Pm. The Pm- 
dependence of rms magnetic field normalised to the equipartition 
field strength, listed in Table [T] shows a much weaker trend. This 
is to be expected as Bcq is proportional to the rms velocity which, 
on the other hand, is a produced by the magnetic field itself. Since 
the parameter range of our simulations is rather limited, no defi- 
nite conclusions can be drawn. However, taking the results at face 
value, it appears that Pm, not Cm, is the parameter that determines 
the saturation le vel in the periodic zero net flux case. Recently, 
iFromand ilOld) reached the same conclusion independently for 
the case of Pm = 4. According to our results, the catastrophic 
quenching and the diffusivity-limited MRI length scale scenarios 
would be ruled out. Although there is the possibility that our calcu- 
lations have not been run long enough, the results seem to suggest 
the small-scale dynamo being harder to excite as Pm decreases. 



3.1.2 Vertical field case 

We find that the saturation behaviour is markedly different when 
vertical field boundary conditions are applied (Table[T]and Fig.|3]l. 
The saturation level of turbulence depends only weakly on the mag- 
netic Prandtl number: the difference of the values of itrms between 
Pm = 0.1 and Pm — 10 cases is roughly 50 per cent. Further- 
more, the Mach number decreases as function of Pm, the trend be- 
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Figure 4. Horizontally averaged horizontal magnetic fields Bx (top panel) 
and By (middle) for Run A3 with Cm = 3 ■ 10^ and Pm = 5. The lower 
panel shows the square of the rms-value of the total magnetic field. 

ing weaker but opposite to the periodic case. This is Hkely caused 
by the increase of viscosity by two orders of magnitude rather than 
the intrinsic dependence of the MRI on Pm. This conjecture is sup- 
ported by the saturation values of the magnetic fields which are 
independent of Pm (lower panel of Fig. [3j- The runs in Set B, 
however, seem to fall into two distinct regimes of magnetic field 
strength, where the magnetic energy differs by roughly a factor 
of two. The reason for this apparent discrepancy is that a differ- 
ent mode of the large-scale magnetic field is excited in the dif- 
ferent branches (see below). Similar behaviour of the large-scale 
dynamo has previously been s een in isotropically forced turbulece 
( [Brandenburg & Dobledl 20021) . 

3.2 Large-scale magnetic fields 

In the runs with periodic boundaries we occasionally see the 
emergence of large-scale magnetic fields with a sinusoidal depen- 
dence on z (see F ig. [4), i.e. fc/fci = 1, in accordance with ear- 
lier investigations jLesur & Ogilvieir2008h . Similar large-scale dy- 
namos have recently been reported from nonhel ically forced turbu- 
lence with shear where^the MRI is absent (e.g. lYousef et al ] |2008l : 
[Brandenburg et alj2008h . As in the forced turbulence case a strong 
large-scale field is not present at all times and the fields undergo 
apparentl y random sign changes that are not fully underst ood (see, 
however. Lesur & Ogilvid[2008l : [Bra ndenburg et al."2008'). The in- 
termittent nature of the large-scale fields could also explain the ap- 
parent lack of catastrophical quenching of the time averaged mean 
magnetic field (see TablelTJ. 

In the vertical field runs a strong large-scale dynamo is al- 
ways excited continuously. The two branches of solutions that are 
visible in the total magnetic energy (Fig. [3} are due to different 
modes of the large-scale field. This is illustrated in Fig. [5] where 
the horizontally averaged horizontal magnetic field components are 
shown as functions of time for Run B 1 . As is common for dynamos 




so 40 80 aO 100 




30 40 60 BO 100 




so 40 60 60 100 



Figure 5. Same as Fig.[4|but for Run Bl with Cm = 1.5 ■ 10^ and Pm = 
10. 



with strong shear, the streamwise component of the magnetic field 
is much stronger than the cross-stream one. Although the initial 
condition of the magnetic field is the same in all runs, the large- 
scale field which develops in the non-linear stage can choose any 
of the available wavenumbers consistent with the vertical bound- 
ary condition Bx — By = 0. In practice, the dominant large-scale 
component is fc/fci = 1 or k/ki = ^ in our simulations. The 
large-scale dynamo tends to accumulate energy at the smallest pos- 
sible wavenumber ( Brandenburg 2001), i.e. the largest spatial scale. 
However, if the dominant mode is on som e intermediate scale ini- 
tially , those modes can also be long-lived ^Brandenburg & DobleJ 
[2OO2I) . Ultimately the large-scale field evolves towards final satura- 
tion where the largest possible scale dominates which was seen in 
IjBrandenburg & Dobler (2002) and in some of our runs (cf. Fig.O. 
The fact that the magnetic energy in Runs B2, B6, B7, and B8 is 
smaller is due to the fact that the large-scale field is predominantly 
of the k/ki = 1 flavour, and that final saturation of the large- 
scale magnetic field has not yet occured. ILesur & Ogilvid ( [200 Sh 
found that the toroidal large-scale magnetic field generated in their 
simulations is close to that yielding the maximum growth rate for 
an m = 1 non-axisymmetric instability. Using their notation we 
find a similar result so that Byky/ {—Sy^ 5/ 12^ nop) ~ 0.6 for 
k/ki = 1 and 1.2 for k/ki = ^, using ky/ki = 1 for the m = 1 
mode. However, the full signifigance of this result is as of yet un- 
clear. 

Although the source of the turbulence and the nature of the 
dynamos (kinematic vs. nonlinear) is differe nt between the non- 
helically forced turbule nce simulations (e.g. [Yousef et al.[ [20081 : 
[Brandenburg et al1[2008b and the non-stratified MRI runs such as 
those presented here, it is conceivable that the large-scale field gen- 
eration mechanism is the same. Since the periodic system is homo- 
geneous, the cause of the large-s cale fields cannot be the g-effect 
of m ean-field dynamo theory ( [Moffatj [ 19781 : [Krause & Radled 
[l980l) . which is in simple systems proportional to the density 
gradient or the turbulence inhomogeneity due to boundaries (e.g. 
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Figure 6. Horizontally averaged kinetic helicity 'H from Run B7. The inset 
shows the volume averaged rms-value of H. The shaded area denotes the 
error estimates. 



Figure 7. Viscosity parameter ass a function of Pm for the runs listed 
in Table[T] The dotted lines show ass = const = 6 ■ 10~^ and QgS <^ 
Pjjj2.o j^jjj. reference. 



iGiesecke et all l2005l : iKapvla et alj l2010iJ) . However, a fluctuat- 
ing a with zero me an can also drive a la rge-scale dynamo when 
shear is present (e.g. Vishniac & Brandenbu rg 1997; Sokolov 199_2l 
ISilant'evll2000l : IProctorl 120071) 7 This is the most likely source of 
the large-scale magnetic fields in the present case. Furthermore, it 
is possible that the shear-current and Q, x J-effects can drive a 
large- scale dynamo iRadler 1969l: lRogachevskii & K leeorin 20031 
l2004h . altho ugh present evidence from numerical models does not 
support this jSrandenburg et alj2008l) . 

In the VF runs the impenetrable stress-free ^-boundaries make 
the turbulence inhomogeneous near the boundary. This leads to 
the generation of mean kinetic helicity Hlz) = oj • u, where 
U3 = V X u is the vorticity. The quantity H is important, be- 
cause t he mean-field g-effect i s, in simple settings, proportional to 
it (e.g. iKrause &" RadleJll980l) . Such contributions, however, will 
not show up in volume averages because the sign of the helicity, and 
thus of the Q-effect, are different near the different boundaries. Fig- 
ure[6]shows the horizontally averaged kinetic helicity for Run B7. 
Here we average also in time over the saturated state of the run. 
In most of the volume the kinetic helicity is consistent with zero, 
although there are regions close to the boundaries where non-zero 
mean values are present. The rms-value of T-i, however, is at least 
five times greater than its mean (see the inset of Fig.[6l(. Note also 
that the normalization factor contains the integral scale k\ . A more 
proper definition would be to use the wavenumber where turbulent 
energy peaks which is likely at least a factor of few greater than fci . 
Thus our estimates for the normalised helicity can be considered as 
upper limits. The rather small values of mean helicity and the dom- 
inance of fluctuations suggest that the generation mechanism of the 
large-scale fields could indeed be the incoherent a~shear dynamo. 
However, a conclusive answer can only be obtained by extracting 
the turbulent transport coefficients an d by perform ing mean-field 
modeling of the same system (see e.g. lGresseill201(]|) . 



terandenburg et alj2004h 



where 



{UxUy) = (UxUy) — {UxUy), 



is the Reynolds stress and 

M^y = (b^by) = (B^By) ^ (B ^B y) , 



(16) 



(17) 



(18) 



the Maxwell stress, and where the angular brackets denote volume 
averaging. Here we decompose the velocity and magnetic field into 
their mean (U, B), taken here as the horizontal average, and fluc- 
tuating (u, b) parts. The mean velocities show no systematic large- 
scale pattern and the remaining signal U ~ ©(O.OSurms) is likely 
a residual of averaging over a finite number of cells. The contri- 
bution of mean flows to the angular momentum transport and the 
dynamo process is thus likely to be negligible. 

For the runs in Set A we find essentially the same scaling, 
consistent with Pm'^", with magnetic Prandtl number as in the 
case of the turbulent kinetic and magnetic energies, see Fig. [7] 
This is consistent with the mixing length estimate of turbulent 
viscosity which is p roportional to the turbulence intensity (e.g. 
ISnellman et alj|2009l) . The numerical values of Qgs decrease from 
« 10"^ for Pm = 5, to ass « 10"'' for Pm = 1. In Set B, on 
the other hand, ass is essentially independent of magnetic Prandtl 
number. The value of ass is consistently of the order of 6 ■ 10""^, 
which is significantly greater than that found in mns with periodic 
boundaries. Here the qualitative behaviour of ass resembles that 
of the turbulent kinetic energy, whereas the two different dynamo 
modes seen in magnetic energy are not visible in the angular mo- 
mentum transport. 



3.3 Angular momentum transport 

The main effect of turbulence in astrophysical disks is to enhance 
diffusion which enables efficient accretion. In accretion disk theory 
it is customary to parametrise the turbulent viscosity Vt in terms 
of the Shakura-Sunyaev vi scosity parameter ass, wh ich relates ut 
with the local gas pressure jShakura & Sunvaevlll973h . 

We define the Shakura-Sunyaev viscosity parameter as 



3.4 Discussion 

A possible clue to understanding the convergence problem in 
zero net flux simulations comes from MRI models with den- 
sity stratification: in them the level of turbul ence does converg e 
when the Reynolds numbers are increased jPavis et al.l |2010[) . 
even with perfect conductor of periodic bound aries. Furthermore, 
such setups exhibit a large-scale dynamo (e.g. [Brandenburg et al] 
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ll995l : IStone et al.lll996l : lGressell2oToh wher e the magnetic helicity 
changes sign at the midplane jGressellbOlOl) . 

Recent numerical results from a different setting suggest tha t 
a diffusive flux of magnetic helicity also exists jMitra et a l. 2010). 
Such a flux can alleviate catastrophic quenching by transporting 
oppositely signed magnetic helicity to the midplane where anni- 
hilation occurs. This could explain the successful convergence of 
the stratified MRI runs. In the non-stratified case with periodic or 
perfectly conducting boundaries, however, no net flux of magnetic 
helicity occurs and the large-scale dynamo can be catastrophically 
quenched, shutting off the MRI. When a flux is allowed by chang- 
ing to vertical field boundary conditions, this limitation is removed 
and the large-scale dynamo can operate without hindrance. How- 
ever, this hypothesis requires further study and more careful analy- 
sis of the helicity fluxes that we postpone to a future publication. 



4 CONCLUSIONS 

We present three-dimensional numerical simulations of the mag- 
netorotational instability in an isothermal non-stratified setup with 
zero net flux initially. Using fully periodic boundaries, that do not 
allow the generation of a mean toroidal flux or magnetic helicity 
fluxes out of the syst em, we encounter the convergence problem 
jpromang et alj[2007l) of the MRI: turbulent kinetic and magnetic 
energies, and the angular momentum transport increase approxi- 
mately proportional to the magnetic Prandtl number. Intermittent 
large-scale magnetic fields are observed in the periodic runs. In- 
creasing the Reynolds numbers moderately at a given Pm does not 
appear to markedly change the results in the saturated state. 

When vertical field boundary conditions, allowing the genera- 
tion of a mean flux and a magnetic helicity flux, are used, the MRI is 
excited at least in the range 0.1 < Pm < 10 for our standard value 
of Cm — 1.5 ■ 10*. We find that the saturation level of the turbu- 
lence and the angular momentum transport are only weakly depen- 
dent on the Prandtl number and that strong large-scale fields are 
generated in all cases. The Shakura-Sunyaev viscosity parameter 
has consistently a value of ass ~ 6 • 10~^ in the vertical field case. 
Exploring even lower values of Pm is infeasible at the moment due 
to prohibitive computational requirements but there are no com- 
pelli ng arguments again st a large-scale dynamo operating at low 
Pm terandenburd 1 200^ . We conjecture that the operation of the 
MRI at low Pm is due to the efficient large-scale dynamo in the sys- 
tem. It is conceivable that the dy namo only work s if magnetic helic- 
ity is allowed to escape (see also lVishniacl200^ or annihilate at th e 
disk midplane due to an internal diffusive flux jMitra et^l201(lh . 
However, measuring the magnetic helicity fluxes in the presence of 
boundaries is difficult due to the fact that they are in general ga uge 
depen dent (e.g. Brandenburg et al. 2002; Hubbard & Brandenburg 
|20l3). 

The current results highlight the close connection between dy- 
namo theory an d the theory of magnetised accretion disks (see also 
iBlackmai] ^20 1 ) and the importance of studying the results in a 
common framework (e.g. Gressel 2010). Clearly, a more thorough 
study is needed in order to substantiate the possible role of mag- 
netic helicity fluxes for the excitation and saturation of the MRI. 
We plan to address these issues in future publications. 



ACKNOWLEDGMENTS 

The authors acknowledge Axel Brandenburg for his helpful com- 
ments on the manuscript. The numerical simulations were per- 
formed with the supercomputers hosted by CSC - IT Center for 
Science in Espoo, Finland, who are administered by the Finnish 
Ministry of Education. Financial support from the Academy of Fin- 
land grant Nos. 121431 (PJK) and 112020 (MIK), are acknowl- 
edged. The authors acknowledge the hospitality of NORDITA dur- 
ing the program "Solar and Stellar Dynamos and Cycles". 



REFERENCES 

Balbus, S. A., Hawley, I.R 1991, ApJ, 376, 214 

Blackman, E.G., & Brandenburg, A. 2002, ApJ, 579, 359 

Blackman, E. 2010, Astron. Nachr., 331, 101 

Brandenburg, A., Nordlund, A., Stein, R. F, & Torkelsson, U. 

1995, ApJ, 446, 741 
Brandenburg, A. 2001, ApJ, 550, 824 

Brandenburg, A. & Dobler, W. 2002, Comp. Phys. Comm., 147, 
471 

Brandenburg, A., Dobler, W., Subramanian, K. 2002, Astron. 
Nachr., 323, 99 

Brandenburg, A., Dintrans, B., & Haugen, N.E.L. 2004, AIPC, 
733, 122 

Brandenburg, A. & Subramanian, K. 2005, Phys. Rep., 417, 1 
Brandenburg, A., Riidler, K.-H., Rheinhardt, M. & Kapyla, P.J. 

2008, ApJ, 676, 740 
Brandenburg, A. 2009, ApJ, 697, 1206 
Cattaneo, R & Hughes, D. W. 1996, Phys. Rev. E, 54, 4532 
Chandrasekhar, S. I960, Proc. Nat. Acad. Sci., 46, 253 
Davis, S. W., Stone, J. M., Pessah, M. E. 2010, ApJ, 713, 52 
Fromang, S. & Papaloizou, J. 2007, A&A, 476, 1113 
Fromang, S., Papaloizou, J., Lesur, G. & Heinemann, T. 2007, 

A&A, 476, 1123 
Fromang, S. 2010, A&A, 514, L5 

Giesecke, A., Ziegler, U. & Rudiger, G. 2005, Phys. Earth Planet. 

Int., 152, 90 
Gressel, O., 2010, MNRAS, 405, 41 

Hawley, J.R, Gammie, C. F & Balbus, S. A. 1995, ApJ, 440, 742 
Hubbard, A., Brandenburg, A., 2011, ApJ, 727, 11 
Korpi, M. J., Kapyla, P. J., Vaisala, M. S. 2010, Astron. Nachr., 
331, 34 

Kapyla, P J., Korpi, M. J. & Brandenburg, A. 2008, A&A, 491, 
353 

Kapyla, P J., Korpi, M. J. & Brandenburg, A. 2010a, MNRAS, 
402, 1458 

Kapyla, P J., Korpi, M. J. & Brandenburg, A. 2010b, A&A, 518, 
A22 

Krause F, Radler K.-H. 1980, Mean-field magnetohydrodynamics 
and dynamo theory (Pergamon Press, Oxford) 

Lesur, G. & Ogilvie, G. 2008, A&A, 488, 451 

Liljestrom, A. J., Korpi, M. J., Kapyla, P. J., Brandenburg, A., 
Lyra, W. 2009, Astron. Nachr., 330, 91 

Mitra, D., Candelaresi, S., Chatterjee, P., Tavakol, R. & Branden- 
burg, A.. 2010, Astron. Nachr., 331, 130 

Moffatt, H. K. 1978, Magnetic field generation in electrically con- 
ducting fluids (Cambridge Univ. Press, Cambridge) 

Proctor M. R. E. 2007, MNRAS, 382, L39 

Radler, K.-H., 1969, Monatsber. Dtsch. Akad. Wiss. Berlin, 11, 
194 



© 0000 RAS, MNRAS 000, 000-000 



8 P. J. Kapyla & M. J. Korpi 



Rogachevskii I., Kleeorin N. 2003, Phys. Rev. E, 68, 036301 
Rogachevskii I., Kleeorin N. 2004, Phys. Rev. E, 70, 046310 
Schekochihin, A.A., Iskakov, A. B., Cowley, S. C, McWilliams, 
J. C, Proctor, M. R. E. & Yousef, T. A, 2007, NJP, 9, 300 
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
Silant'ev N. A. 2000, A&A, 364, 339 

Snellman, J. E., Kapyla, P. J., Korpi, M. J., & Liljestrom, A. J., 

2009, A&A, 505, 955 
Sokolov D. D. 1997, Astron. Reports, 41, 68 
Stone, J. M., Hawley, J. R, Gammie, C. R & Balbus, S. A. 1996, 

ApJ, 463, 656 
Vainshtein, S. I. & Cattaneo, R 1992, ApJ, 393, 165 
Velikhov, E.R 1959, Sov. Phys. JETP, 36, 1398 
Vishniac E. T., Brandenburg A. 1997, ApJ, 475, 263 
Vishniac E. T., Cho J. 2001, ApJ, 550, 752 
Vishniac E. T. 2009, ApJ, 696, 1021 
Wisdom, J., Tremaine, S. 1988, AJ, 95, 925 
Yousef T. A., Heinemann T., Schekochihin A. A., Kleeorin N., 

Rogachevskii 1., Iskakov A. B., Cowley S. C, McWilliams J. C. 

2008, Phys. Rev. Lett., 100, 184501 



